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Abstract 

Dilute multi a cluster condensed states with spherical and axially deformed shapes are 
studied with the Gross-Pitaevskii equation and Hill- Wheeler equation, where the a cluster 
is treated as a structureless boson. Applications to self-conjugate 4A^ nuclei show that the 
dilute Na states of ^■^C to ^'^Ca with J'^ = 0"*" appear in the energy region from threshold 
up to about 20 MeV, and the critical number of a bosons that the dilute Na system can 
sustain as a self-bound nucleus is estimated roughly to be N^r ~ 10. We discuss the charac- 
teristics of the dilute Na states with emphasis on the dependence of their energies and rms radii. 

PACS numbers: 21.10.Dr, 21.10.Gv, 03.75.Hh 
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I. INTRODUCTION 



The molecular-like picture as well as the sinerle-Darticle picture are fundamental view- 
points to understand the structure of light nuclei. It is well known that the 
structure of many states in light nuclei is described successfully by the microscopic cluster 
model, where a group of nucleons is assumed to form a localized substructure (cluster), 
interacting with other clusters and/or nucleons in the nucleus. The predominant cluster is 
the a nucleus, which plays an important role in the cluster model, because it is the lightest 
and also smallest shell-closed nucleus with a binding energy as large as 28 MeV, reflecting 
the strong four-nucleon correlation. Molecular-like states in nuclei are expected to appear 
around the threshold energy of breakup into constituent clusters ^ , because the intercluster 
binding is weak in the cluster states. For example, the ground state of ^Be and the second 
O"*" state of ^^C are known to have loosely bound 2a and 3a structures, respectively, which 
appear around the 2a and 3a thresholds. 

Special attention has been paid to such four-nucleon correlations corresponding to an a- 
type condensation in symmetric nuclear matter, similar to the Bose-Einstein condensation 
for finite number of dilute bosonic atoms such as ^''Rb or ^^Na at very low temperature 
Several authors have discussed the possibility of a-particle condensation in low-density nu- 
clear matter 0,0]. They found that such a condensation can occur in the low-density region 
below a fifth of the saturation density, although the ordinary pairing correlation can prevail 
at higher density. The result indicates that a condensate states in finite nuclear system may 
exist in excited states of dilute density composed of weakly interacting gas of a particles. 
Thus, it is an interesting subject to study the structure of light nuclei from the viewpoint 
of a-particle condensation. 

Recently, a new a-cluster wave function was proposed which is of the A^a-particle con- 
densate type [l^ : 

\^Na) = (C^fW^c), (1) 

(r-1 ■ ■ ■ r^|$^„) oc A |e-K^?+- +^^)0(aO ■ ■ ■ , (2) 

where is the a-particle creation operator, (j){a) denotes the internal wave function of 
the a cluster, is the center-of-mass coordinate of the i-th a cluster, and A presents the 
antisymmetrization operator among the constituent nucleons. The important characteristic 
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in Eq. (jj)) is that the center-of-mass motion of each a cluster is of 5'-wave type with the 
independent size parameter v. Apphcations of the condensate-type wave function to ^^C 
and ^^O showed that the second 0+ state of (E^=7.65 MeV) and fifth 0+ states of ^^O 
(-E^.=14.0 MeV), located around the 3a;- and 4a;-particle thresholds, respectively, are spec- 
ified by the condensate state, which is quite similar to the Bose-Einstein condensation 
of bosonic atoms in magnetic traps where all atoms occupy the lowest S'-orbit.^^ The cal- 
culated root-mean-square (rms) radius for those condensate states is about 4 fm, which is 
much much larger than that for the ground state (about 2.7 fm). As for ^Be, the a-particle 
wave function in Eqs. and taking into account the axially symmetric deformation, 
was applied to investigate the rotational structure of the ground-band states (0+-2+-4"'") 
It was found that the rotational character of the ^Be ground state band is reproduced 



nicely by this wave function with deformation [11 1. 

The above-mentioned theoretical results for ^Be, ^^C and ^^O lead us to conjecture that 
such dilute a-cluster states near the Na threshold may also occur in other heavier 4A^ self- 
conjugate nuclei. The Coulomb potential barrier should play an important role to confine 
such dilute A^a-particle states, as inferred from the analyses of the ^Be, ^^C and ^^O nuclei. 
In fact, the ground state of ^Be, which appears at = 92 keV referring to the 2a threshold, 
exists as a resonant state with very narrow width, due to the Coulomb potential barrier, 
whose height is estimated to be 1 ~ 2 MeV. This self-trapping of a-particles by the Coulomb 
barriers is in contrast to the case of the dilute neutral atomic condensate states, where the 
ato,n.a.e capped by the extent. ,na,„etic«e>d0. 

It is an intriguing problem and of interest to study how many a particles can be bound in 
the dilute nuclear system. Increasing the number of a clusters, the rms radius of the system 
should become gradually larger, because the dilute character must be retained. Then, the 
total kinetic energy of the Na system becomes significantly smaller in comparison with the 
potential energy, similar to the case of bosonic atoms in the condensed state. In addition, 
the height of the Coulomb potential barrier will become steadily lower. This is due to the 
following reasons. The a-a nuclear potential is short-ranged, while the Coulomb potential 
is long-ranged. Increasing the number of the a particles, the repulsive contribution from 
the Coulomb potentials prevails gradually over the attractive one from the a-a nuclear 
potentials, because the gas-like Na system is expanding in such a way that the rms radius 
between two a particles is gradually getting larger. This means that the height of the 
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potential barrier confining the gas-like Na particles becomes lower with increasing N ^ and 
then, there should exist a critical number (A^ = Ncr) beyond which the Na system can 
not sustain itself as a bound nuclear state. Thus, it is interesting to estimate the number 
of Ncr as well as to study the structure of the dilute Na system up to = N^r- Since 
the application of the condensate-type wave function in Eqs. (jT)) and (j2I) to the general Na 
system is not easy and the computational limitation is ~ 6 at most, we need to study 
them with phenomenological models. 

In this paper, the gas-like Na cluster states are studied with the following two approaches: 
the Gross-Pitaevskii-equation approach and Hill- Wheeler-equation approach. The Gross- 
Pitaevskii equation, which is of the nonlinear-Schrodinger-type, was proposed about fifty 
years ago to describe the single-boson motion in the dilute atomic condensate state |l2j |. 
Recent experiments starting from the middle of 1990's have succeeded in realizing such a 
condensate state consisting of 10^-10^ neutral atoms trapped by the magnetic field at very 
low temperature. Many characteristic aspects of the dilute states are described successfully 
by the Gross-Pitaevskii equation Thus, its application to the dilute Na nuclear systems 
is also very promising, and we expect from such a study useful information on the condensed 
states. On the other hand, we here propose also a different approach. Although the Gross- 
Pitaevskii equation is simple and interesting for the study of the structure of the dilute 
A^a system, the center-of-mass motion in the A^a system is not completely removed in 
this framework. The effect of the center-of-mass motion should be non-negligible to the 
total energy and rms radius etc. for small numbers of a bosons. We, thus, formulate in 
this paper the framework describing the dilute A^a systems free from the center-of-mass 
motion, with the approach of the Hill- Wheeler equation. The spherical A^a systems as well 
as the deformed ones, with the axial symmetry, are discussed with this equation. The two 
different approaches, the Gross-Pitaevskii-equation approach and the Hill- Wheeler-equation 
approach, are complementary to one another, and we will obtain useful understanding of 
the dilute multi a-particle states. Starting from ^Be, the structure of the gas-like A^a boson 
systems with J'^ = 0+ is investigated and the critical number Ncr is estimated with those 
two different approaches. 

The paper is organized as follows. We formulate in Sec. II the Gross-Pitaevskii-equation 
approach and the Hill- Wheeler-equation approach for the dilute A^a system. In Sec. HI the 
calculated results are presented, and the characteristics of the dilute A^a states are discussed 
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with emphasize on the N dependence of their energies and rms radii. A summary finally is 
given in Sec. IV. 



II. FORMULATION 

In this section, we formulate the two approaches to study the structure of the dilute 
Na nuclear systems: the Gross-Pitaevskii-equation approach and Hill- Wheeler-equation 
approach. 

A. Gross-Pitaevskii equation for dilute Na nuclear systems 

In the mean field approach, the total wave function of the condensate A^'a-boson system 
is represented as 

N 

ma)^l[cp{ri), (3) 

i=l 

where (p and are the normalized single-o; wave function and coordinate of the i-th a boson. 
Then, the equation of motion for the a boson, called as the Gross-Pitavskii equation, is of 
non-linear Schrodinger-type, 

U{r) = {N-1) J dr' \ip{r')\\2{r' , r) 

+ l(7V-i)(7V-2) I dr"dr'Mr")\'Mr')\'vs{r",r',r), (5) 

where m stands for the mass of the a particle, U is the mean-field potential of ct-particles, 
and V2 (fa) denotes the 2a (3a) interaction. The center-of-mass kinetic energy correction, 
1 — is taken into account together with the finite number corrections. A" — 1, etc. In the 
present study, only the S-wave state is solved self-consistently with the iterative method. 
The total energy of the Na system E(Na) is presented as 



E{Na) = N 



{t) + 1{N- 1){V2) + -1){N- 2){v,) 



(6) 



(^)=(l-^)xMr)|-^V^|^(r)), (7) 
{V2) = {cp{r)ip{r')\v2{r,r')\ip{r)ip{r')), (8) 
(vs) = (^(r)^(0¥^(r")|f3(r-,r',r")|<^(r)¥.(0¥'(r")), (9) 



and the eigen energy e in Eq. Q is given as 

e={t) + {N- 1){V2) + \{N - 1){N - 2){vs). (10) 

The nuclear rms radius in the Na state is evaluated as, 

= ^K)gp + l.7l\ (11) 
{rl)GP={l-j^){v\r'\'p), (12) 

where we take into account the finite size effect of the a particle and the correction of the 
center-of-mass motion. 



B. Hill- Wheeler equation for dilute A^q states 

The framework of the Gross-Pitaevskii approach is simple and useful for the study of the 
structure of the dilute Na systems. The center-of-mass motion, however, is not removed 
exactly in this approach. The effect on the total energy and rms radius etc. is not negli- 
gible for the small-number a-boson systems. We formulate here the Hill- Wheeler-equation 
approach for the axially symmetric Na systems as well as for the spherical ones, in which 
the center-of-mass motion is completely eliminated. Only the Hill- Wheeler equation for the 
deformed Na systems with the axial symmetry is provided in this section, because the limit 
of the spherical systems is included in the deformed case. 

The single a-particle wave function should be mainly in the lowest S state in the gas-like 
Na boson system. Thus, the model wave function of the axially symmetric Na system 
should be described in terms of a superposition of a Gaussian basis with axially symmetric 
deformation (taking the z-axis as the symmetric axis) as follows, 

$M)(iVa) = J2 fi'^u ^3)$^*"*n^i, i^s), (13) 

$(^"*)(z/i, z/3) = / rfi?e™$^^"^*(i?cm; 1^1, i^sMi^i, 1^3), (14) 

N 

= n</'(^-'^l''^3) = $(*"*n^l,i^3)$^^"^n^cm;i^l,^^3), (15) 
i=l 

(t>{r] ui, z/3) = (^^^ exp (-z^ix^ - uiy^ - v^z^^ , (16) 

^(^"^n^cm; yi, yz) = <l>{Rcm; Nui, Nus), (17) 
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where (Rem) denotes the coordinate of the z-th a boson (the center-of-mass coordinate of 
the Na system), and z/i (z/s) presents the Gaussian size parameter of the x and y directions 
{z). The wave function in Eq. (|13p is totally symmetric under any exchange of two a bosons. 
The angular-momentum-projected total wave function for dilute Na systems free from the 
center-of-mass motion is given as 

^^}M\Na) = E /X^^i, ^^3)$SVi, ^3), (18) 



$S)(z/i,z/3) = / dcos9di,{9)ny{9)<^^''''\u^,U3), (19) 

where TZy denotes the rotation operator around the y axis. Setting z/3 = ui in Eq. (fTHjl . we 
obtain the wave function for the spherical Na system, where only = 0"*" state is allowed. 
The matrix element of a translational invariant scalar operator O with respect to the angular 
momentum projected Na wave function in Eq. (fT^ is evaluated as 

(4mLo(^1,^3)|0|$SLo(^i',^3')) 

= I dcos<o(^)('^'^*"*H^l,^3)|7^,(^)C)|$('"*)(^l^^3')), (20) 

.cos<o(^)7,,r^7#^^^^ (21) 

The total Hamiltonian of the A^a-boson system is given as 

N 

'H = ^U-Tcm + ^V2{ri,rj) + ^ V3{ri,rj,rk), (22) 

1=1 i<j i<j<k 

where V2 and V3 denotes, respectively, the 2a and 3a interactions. The kinetic energy of the 
center-of-mass motion [Tcm) is subtracted from the Hamiltonian. 

The equation of motion for the dilute A^a-boson states are given in terms of the Hill- 

n 

Wheeler equation jl3| : 

Y: {($rVl,^3)|H-i?|$f*^(^l',^3'))}^(^l',^3') =0, (23) 

where Ti is the total Hamiltonian of the A^a-boson system in Eq. ()22j) . The coefficients fj 
and eigen energies E are obtained by solving the Hill- Wheeler equation. In the numerical 
calculation, the Gaussian size parameters i>i and z/3 are discretized and chosen to be of 
geometric progression, 

u['^ = {l/b?)\ 6f) = 6«rr, fc = l~A:_, (24) 
-f^ = (lM"T' bf^ = b?rf-\ ir=l^ir_., (25) 



The above choice of the Gaussian range parameters is found to be suitable for describing 
the dilute Na states. 

The nuclear rms radius measured from the center-of-mass coordinate in the Na state, is 
expressed as, 

= /(^iW + m^, (26) 

{rl)HW - {^T\Na)\^Y.ir^ - Rcmn^THNa)), (27) 

where we take into account the finite size effect of the a particle. The rms distance between 
two a particles is given as 

Thus, it is proportional to the rms radius of an ct-particle from the center-of-mass coordinate. 



C. Effective a-a potentials. 

In the present paper, we use two kinds of effective potentials: the density-dependent po- 
tential and phenomenological 2a plus 3a potential. They are applied to the Gross-Pitaevskii 
equation and the Hill- Wheeler equation for the study of the structure of dilute Na states. 



1. Density- dependent potential 

The density dependent potential consists of the Gaussian-type a-a potential including 
a density-dependent term, which is of similar form as the Gogny potential (known as an 
effective NN potential) used in nuclear mean-field calculations. 



V2{r,r') = -uoexp [-0.7^(r - r')^\ - 130 exp [-0.475^(r - r'f 
+{4Try gS{r - r')p — — + vcoui{.r, r'), 



(29) 



where the units of V2 and r arc McV and fm, respectively, and p denotes the density of the 
Na system. The folded Coulomb potential vcoui is presented as 



vcoui{r,r') 



4e2 



\r — r 



— erf(a|r — r'|). 



(30) 



8 



The Gaussian-potential part in Eq. ()29|) is based on the Ah-Bodmer potential jlj|, which 
is known to reproduce well the elastic a-a scattering phase shift up to about 60 MeV for 
= 500 MeV. On the contrary, here, the two parameters, Vq and g, are chosen so as to 
reproduce well the ex per imental energy (£^3^^^=0.38 MeV) and the calculated rms radius (4.29 
fm by Tohsaki et al. [lO|) for the 02^ state of ^^C by solving the Gross-Pitaevskii equation in 
Eq. (0)). The results are t;o=271 MeV and g=1650 MeV-fm^, where the calculated energy and 
rms radius for the 3a system are 0.38 MeV and 4.14 fm, respectively, which are discussed in 
Sec. III. Although the choice of the phenomenological potential in Eq. ()29|) is rather rough, it 
is interesting to study systematically the structure of the Na condensate states as a function 
of N. 



2. Phenomenological 2a and 3a potential 

There are many phenomenological 2a potentials proposed so far. Since the a particle 
is treated as a point-like boson in the present study, we will use a 2a potential taking 



3, 



into account the Pauli blocking effect. The typical potential is the Ali-Bodmer one 
which is used frequently in the structure calculation; V2{r,r') = 500 exp [— 0.7^(r — r')^] — 
130 exp [— 0.475^(r — r')^] + Vcouiix, f'), where Vcoui denotes the folded Coulomb potential 
given in Eq. (j30p . The strong repulsion in the inner region prevents the 2a particles from 
approaching one another. It is, however, found that the potential is not suitable to describe 
the property of the compact shell- model-like structure of the ^^C ground state with the 
3a; boson model However, we may use it for the dilute A^a states. The Ali-Bodmer 
potential, however, has the following three unfavorable properties for the present calculation. 

First is that the potential does not give the experimental resonant energy of the ^Be 
ground state (i?2a=68 keV vs. £'2^^=92 keV ), although the a-a scattering phase shift is 
reproduced nicely up to about 60 MeV. The second is that applying the potential to the 
3a boson system the lowest energy state obtained corresponds to a relatively compact 3a- 
structure state, although the condensate state appears around the 3a threshold. According 
to the stochastic variational calculation ^|, for example, the calculated energy and rms 
radius are, respectively, E^'^ = —0.62 MeV and J (r^)=3.15 fm The results indicate 



that the Ali-Bodmer potential is not adequate for describing the dilute 3a-boson state of 
^^C. The third unfavorable point comes from the fact that the strong repulsive character in 
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the short-range region of the Ah-Bodmer potential (~ 400 MeV) leads us to treat exactly the 
short-range correlation between the la bosons. The treatment is difficult and needs time- 
consuming numerical calculations for solving general iVa-boson systems, even if employing 
modern numerical methods for many body systems. 

We construct here an effective la potential with a weak repulsive part (soft core), which 
overcomes the above-mentioned three unfavorable properties. The use of such a soft-core- 
type la potential is suitable in the present study, because we discuss the gas-like Na states, 
where the Pauli blocking effect between two a particles is considerably weakened. The 
important point for the determination of the potential parameters, other than the condition 
of reproducing the experimental resonant energy for the ^Be ground state, is that the a-a 
wave function for the resonant state, should have a loosely-bound structure of the two a 
particles. According to many structure calculations of ^Be, the amplitude of the radial part 
of the a-a relative wave function must be small in the inner region and have a maximum 
value around r = 4 fm (where r is the relative coordinate between the two a particles). This 
condition ensures that the ground state of ^Be has a dilute 2a structure. With a careful 
search of the potential parameters, we determined the effective 2a potential as follows, 

v^{r, r') = 50 exp [-OA'^{r - r')^] - 34.101 exp [-0.3'^{r - r'f] + vcouiir, r'), (31) 

where the units of V2 and r are MeV and fm, respectively, and Vcoui is the folded Coulomb 
potential. The calculated resonant energy of ^Be is i?2a=92 keV, in agreement with the ex- 
perimental data {E2a^=92 keV). Figure ^ shows the radial part of the relative wave function 
between the 2a clusters for the resonant state. We see that the amplitude is relatively small 
for r = ~ 2 fm, and has a maximum value around r=4 fm. In spite of the fact that the 
Ali-Bodmer potential gives a wave function which is almost zero for r = ~ 1 fm, the small 
but finite amplitude around r = ~ 2 fm should hardly give any effect for dilute multi a 
cluster states. Applying this effective potential to the 3a- and 4a-boson systems indicates 
that we get the desired dilute 3a- and 4a-structure states for ^^C and ^^O, respectively, 
within our framework, as shown below. 

In the present study, we introduce the phenomenological 3a potential (v^) with repulsive 
character, as given in Ref. 

vs{r, r', r") = 151.5 exp {-0.15 [(r - r'f + (r' - r")^ + (r" - rf] }, (32) 
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where the units of vs and r are MeV and fm, respectively. This potential has been used in the 
3a and 4a orthogonally condition model (OCM for calculations of ^^C and ^^O so as to 
reproduce the ground-state energies with respect to the 3a and 4a thresholds, respectively 
[l8|. In the model, the Pauli principle is taken into account in the relative wave function 
between two a particles, and a deep attractive potential is used for the a-a potential. Thus, 
the OCM is able to describe not only the shell-model-like compact states but also the dilute 
gas-like states. According to the results, the repulsive 3a potential gives a large effect to 
the ground-state energies of ^^C and ^^O with the compact Na structure, while its effect is 
very small for dilute 3a and 4a states. A non-negligible effect, however, can be expected 
in large-number A^a dilute systems, if we take into account the fact that the contribution 
of the binding energy from the 3a potential, proportional to A^(A^ — 1)(A^ — 2)/6, raises 
strongly with increasing A^. 

The reason of why we introduce the repulsive 3a potential is given as follows. Let us 
define here the total kinetic energy and two-body potential energy in the dilute A^a system 
as (T) and (V2), respectively. In case of the ^Be ground state, the experimental energy is 
E = (T) -|- (V2) ~ 0.1 MeV with respect to the 2a threshold, where (T) ((V2)) is positive 
(negative). For an arbitrary dilute A^a-boson system, the total kinetic energy and two-body 
potential energy are given as (T) ~ A^ — 1 and (V2) ~ A^(A^ — l)/2, respectively, where 
the center-of-mass kinetic energy is subtracted. Increasing the number of the a particles, 
therefore, the potential energy prevails over the kinetic energy, and then, the system falls 
gradually into a collapsed state. This indicates that something like a density- dependent force 
with the repulsive character is needed to avoid the collapse in the large-number A^a-boson 
system. The present repulsive 3a potential in Eq. (jH^ plays a role similar to the density- 
dependent force. On the other hand, the density-dependent potential given in Eq. (j^^ is also 
used when solving the Gross-Pitaevskii equation and the Hill- Wheeler equation. The reason 
of why the density- dependent potential is introduced there is the same as that discussed 
here. 
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III. RESULTS AND DISCUSSION 



A. Application of the Gross-Pitaevskii equation to Na systems 

The Gross-Pitaevskii equation is solved with the two different types of effective a-a 
potentials: 1) the density-dependent potential [see Eq. (j^ ] and 2) the phenomenological 
2a potential with the 3a potential [see Eqs. (jM]) and (1221)] • First of all, we will discuss 
the results with the density-dependent potential and then those with the phenomenological 
potentials. 

The calculated total energies of the Na systems measured from the Na threshold are 
demonstrated in Fig. |21as well as the calculated nuclear rms radii defined in Eq. where 
the density-dependent a-a potential is used. The total energy and the rms radius are getting 
larger with increasing N . This means that the system is expanding steadily with increase 
of N . In comparison with the rms radius of the ground state of each nucleus with the 
empirical formula {\J (r^) = 1.2^4^^^ fm), the results from the Gross-Pitaevskii equation 
are much larger than those for the ground states. Thus, the states obtained here can be 
identified with the dilute Na states. It is noted that they are obtained naturally from 
the Gross-Pitaevskii equation with the density-dependent a-a potential, whose parameters 
were determined so as to reproduce well the experimental energy (i?3^^=0.38 MeV) and the 
calculated rms radius by Tohsaki et al. (4.29 fm) for the state of ^^C. 

In order to study the structure of the dilute Na states, it is instructive to see the single 
a potential defined in Eq. Figures inta)~(j) show the ones for the 3a ~ 11a systems. 
Let us first discuss the 3a and 4a cases. The remarkable characteristics of the potentials 
can be presented as follows: 1) the almost flat behavior of the potential in the inside re- 
gion, and 2) the Coulomb potential barrier in the outer region. The appearance of the flat 
potential region is very impressive, if one takes into account the fact that the two-range 
Gaussian term (attractive), density-dependent term (repulsive) and the Coulomb-potential 
term (repulsive) contribute significantly to the single a potential in the inside region [see 
Figs. Efa) and (b)]. Table IJ shows the calculated single a particle energy and contributions 
from the kinetic energy, two-range Gaussian term in Eq. (j^^ . density-dependent term in 
Eq. ()29|) and the Coulomb potential. It is found that the kinetic energy in the 3a and 4a 
systems is not negligible but small in comparison with the two-range- Gaussian term and/or 
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Coulomb potential energy. This indicates that the Thomas-Fermi approximation, neglecting 
the kinetic energy term, is roughly realized in the system. Consequently, there appears the 
fiat potential region in the single a potential, whose behavior is similar to the dilute atomic 
condensate state trapped by the magnetic fields at very low temperature |7]. On the other 
hand, the appearance of the Coulomb potential barrier plays an important role in confining 
the a bosons in the inside region. It is noted that the barrier comes out naturally from the 
self-consistent calculation of the Gross-Pitaevskii equation. The single a particle energy for 
the 3a and 4a systems are smaller than the Coulomb potential barrier. This means that 
the dilute states are quasi-stable against a decay. 

Increasing the number of a bosons beyond = 5, the following interesting features 
can be seen in the single a particle potentials: 1) the depth of the fiat potential becomes 
shallower, and its range is expanding to the outer region faster than oc A^^/^, and 2) the 
height of the Coulomb potential barrier is getting lower and almost disappears at around 

= 10. The first behavior of the potentials means that the A^a system is infiating with 
increase of A^, as inferred from the behavior of the rms radius shown in Fig. |21 The reason 
of why the depth of the potential becomes shallower is related to the fact that the single 
a particle energy is getting larger. It is given as follows: From Table lU the contribution 
from the kinetic energy to the single a particle energy becomes smaller with increase of A^, 
refiecting the infiation of the dilute A^a system. This fact indicates that the Thomas-Fermi 
approximation is getting better, and then, the single a particle energy is given approximately 
as the contribution from only the potential energies. The potential energies consist of three 
contributions, namely, the two-range Gaussian term, density-dependent term, and Coulomb 
potential [see Eq. (j^^ ]. where the first and second ones are short-range, and the third is 
long-range. When the dilute A^a system is infiating with A^ and the distance between two 
a bosons is getting larger, the contribution from the long-range potential should overcome 
steadily the one from the short-range one. Consequently, the potential depth (single a 
particle energy) becomes gradually shallower (larger) and the height of the Coulomb barrier 
is getting lower, with increase of A^. The increase of the single a particle energy means that 
the total energy of the dilute A^a system becomes larger with A^. The present results are 
consistent with the behavior of the A^ dependence of the total energy (see Fig. 12)). 

The second interesting behavior of the single a particle potential with A^ > 5 is that there 
exists a critical number of a bosons, A"cr, beyond which the system is not confined anymore. 
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as mentioned in Sec. I. Around = 10, the Coulomb barrier has almost disappeared, and 
the single a particle potential is nearly flat up to the outer region, although the most outer 
region, dominated by only the Coulomb potential, is falling to zero (not illustrated in Fig. 01). 
Thus, we can roughly estimate the critical number as Ncr ~ 10, namely, ^°Ca. 

It is interesting to compare the above results with those obtained by using the Gross- 
Pitaevskii approach with the phenomenological 2a and 3a potentials given in Eqs. ()31|) and 
(jH^ . respectively. The calculated total energies and rms radii for A^a systems are shown in 
Fig. 0] We find that they are in good agreement with those of the case with the density- 
dependent potential in the Gross-Pitaevskii equation within about 10% except for the A^ = 3 
and 4 cases. Table ITTl shows the single a particle energies and contributions from the kinetic 
energy and potential energies. Comparing them with those in Table HJ the A^ dependence 
of the kinetic energy. Coulomb potential energy and the sum of other potential energies 
is rather similar to the one with the density-dependent potential in the Gross-Pitaevskii 
equation except for the A^ = 3 and 4 systems. This is surprising if taking into account 
the fact that the forms of the two kinds of the effective a-a potentials are quite different 
from one to another. These results might indicate that the dilute A^a systems with A^ > 5 
do not depend strongly on the details of the effective a-a potential, while those with a 
small-number- A^a (A^ = 3 and 4) are sensitive to the potential. 

The single a particle potentials are shown in Fig. |S| For the 3a and 4a systems, we see 
the almost flat potential behavior in the inside region, and the Coulomb potential barrier 
in the outer region, while increasing the number of a bosons from A^ = 5, the depth of the 
flat potential becomes shallower, and its range is expanding to the outer region, and the 
height of the Coulomb potential barrier is getting lower and almost disappears at around 
A^ = 10. The qualitative potential behaviors are almost the same as those in case of the 
density-dependent potential with the Gross-Pitaevskii equation (see Fig. E)). In fact, the 
behavior of the single a particle potential in Fig. El is in good agreement with that in Fig. |H1 
within about 10 %. On the other hand, we can conjecture the critical number Ncr from the 
behavior of the single a potentials in Fig. which is estimated roughly as Ncr ~ 10, the 
result being the same as that of the case with the density-dependent potential. 
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B. Na systems in the Hill- Wheeler-equation approach 

The Gross-Pitaevskii equation is simple and useful for the study of the structure of the 
dilute Na system, as seen in the previous subsection. The center-of-mass motion in the 
Na system, however, is not completely removed in this framework. The Hill- Wheeler- 
equation approach is free from the center-of-mass motion of the Na system. The effect 
of the center-of-mass motion is expected to be non-negligible for the total energy and rms 
radius etc. in small-number a-boson systems, in particular, = 3 and 4. Thus, it is 
important to study the dilute Na systems with the use of the Hill- Wheeler equation. The 
two approaches should give the same results for dilute Na systems with rather large number 
of A^. In this subsection, we first demonstrate that the Hill- Wheeler-equation approach is 
useful to describe the dilute Na states with the spherical shape, and the results are briefly 
compared with those from the Gross-Pitaevskii equation. Then, the approach is applied to 
the deformed Na systems {J'^ = 0^) with axial symmetry. We use the phenomenological 
2a and 3a potentials of Eqs. pijl and (j32|) as the effective aa potential. 

1. Spherical Na systems 

The calculated total energies and rms radii for spherical Na states are illustrated in 
Fig. ini in which we show those obtained by solving the Gross-Pitaevskii equation with the 
same effective a-a potential for comparison. We found that the calculated total energies in 
the Hill- Wheeler approach are almost the same as those in the Gross-Pitaevskii approach. 
This indicates that the center-of-mass kinetic energy correction in Eq. (0]) is a very good 
approximation to remove the effect in the total energy. The non-negligible deviation between 
the two frameworks, however, can be seen in the rms radii, in particular, we have about 20 % 
deviation for the 3a system, although that is getting smaller with increasing and it is 
almost zero in the 7a ~ 9a systems (the discrepancy in the 10a and 11a systems will be 
discussed later). The results tell us that the center-of-mass correction for the rms radius 
in the Gross-Pitaevskii-equation approach [see Eqs. ()11|) and (fT^ ] is not very good for the 
3a and 4a systems, while it is relatively good for the A^a > 5a systems, and the difference 
between the two approaches is less than about 10 %, and it is diminishing with increasing 
N. 
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In Fig. iniwe see that the deviation of the rms radius appears again in the lOa and 11a 
systems. The reason of why this deviation occurs can be understood, recaUing the following 
facts: According to the results of the Gross-Pitaevskii equation with the phenomenological 
la and 3a potentials, the 10a system is the critical one which can not exist as a nuclear 
state, as mentioned in the previous section. In such a critical system, the behavior of the 
wave functions in the outer (inner) region is very sensitive (not very sensitive) to how to 
solve the equations and to obtain the wavefunctions. The value of the rms radius (total 
energy) is generally sensitive (not very sensitive) to the behavior of the wave function in the 
outer region. Thus, the re-occurrence of the deviation in the rms radii indicates that the 
10a system is critical in the Hill- Wheeler-equation approach. 

The above results show us that the Hill- Wheeler-equation approach is very useful to 
describe the dilute states as well as the Gross-Pitaevskii-equation approach. Thus, we 
can apply the approach to the deformed system with the axial deformation. 



2. Deformed Na systems with = 0^ 

The Na states {J'^ = 0~^) with the axial deformation are obtained by solving the Hill- 
Wheeler equation in Eq. (j^Hj) with the Gaussian size parameter set A given in Table IIIIl 
Figure [7| illustrates the calculated energies, nuclear rms radii defined in Eq. (|11|1 and the rms 
distances between 2a bosons. The detailed values are shown in Table IIVI The rms distance 
between 2a bosons (6 ~ 11 fm) for = 3 ~ 12 is considerably larger than that of the ®Be 
ground state (~ 4 fm). The result indicates that the A^a states obtained here are of very 
dilute A^a structure. 

The total energy of the A^a state increases gradually with A^, although those for A^=3 
and 4 are not changed very much. The latter is in contrast to the results of the spherical 
case together with those for the Gross-Pitaevskii-equation approach (see Figs. El El and Ej). 
In the 3a and 4a systems, the total energies and nuclear rms radii are respectively given as 
follows; E = —0.01 and 0.13 MeV, and (r^)=3.73 and 3.90 fm. The values are in good 
correspondence to those for the condensate a cluster states discussed by Tohsaki et al. |lO| , 
where their calculated results are E™'=0.5 and -0.7 MeV (vs. ^'^^^=0.38 and -0.44 MeV) 
and ^ (r^)=4.29 and 3.97 fm, respectively, for the dilute 3a and 4a states. It is instructive 
here to compare the present results with those in the spherical A^a states in order to estimate 
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the effect of deformation. The comparison of the total energy and rms radius is given in 
Table HVl The energy gains (reductions of rms radius) due to the deformation in the 3a, 4a 
and 5a systems are, respectively, 0.6, 1.5 and 1.7 MeV (11 %, 18 % and 19 %), while this 
is getting smaller for > 6. Thus, we found that the deformation effect is significant for 
relatively small-number a systems, and shows a good correspondence of our results for the 



3a and 4a states with those by Tohsaki et al. [l^. The present result that the dilute Na 
states {J'^ = 0+) may be deformed is natural if taking into account the fact that a gas-like 
Na state with relatively small number can easily be deformed. If it is right, there may exist 
dilute Na nuclear states with = 2^ and 4"'" etc. In fact, a candidate of the dilute 3a 
state w,th . 2. ,s observed at ^^3.3 MeV „,easu..ed .o-n tKe Sa threshold Q. The 
dilute multi a cluster states with non-zero angular momentum will be discussed elsewhere. 

From Fig. [7(a), we notice that the gas-like Na states with > 5 appear above the 
Na threshold, not close to it, in contrast to the fact that the dilute 3a and 4a states are 
located in the vicinity of their respective thresholds. This is a non-trivial result of the 
present calculation. The reason of why the energy of the dilute Na states increases with 

is given as follows. Let us consider the effective energy of the "2a" system in the dilute 
A^a state, which corresponds to the quantity of the total energy of the A^a state divided 
by the number of a pairs which is given in Table IIVI Since the dilute character allows us 
to neglect approximately the kinetic energy, the effective energy of the "2a" system could 
be given mainly as the sum of the 2a nuclear potential energy and its Coulomb potential 
energy (the contribution from the 3a potential may be neglected approximately because of 
the very small amount). Figure |H1 shows the 2a potential used in the present study, where 
we also show the 2a nuclear potential [see Eq. (jHH)] and its Coulomb potential [see Eq. (jHUj) ]. 
The attraction of the former is largest around 4 fm and is negligible beyond about 7 fm, 
while the Coulomb potential is long-ranged and its repulsion is substantial even beyond 7 
fm. Increasing the number of the a particles, the gas-like A^a system is expanding and the 
rms distance between two a particles is also becoming larger (6 ~ 10 fm) in the case of 
A^ > 5, as shown in Fig. [3 Then, the attractive contribution from the 2a nuclear potential 
to the effective "2a" energy should become noticeably smaller with increasing A^, reflecting 
the short-range (attractive) behavior of the potential. The repulsive contribution from the 
Coulomb potential, however, should not be much smaller than in the case of the 2a nuclear 
potential, and dominate over the contribution from the nuclear potential, because of the 
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long-range (repulsive) behavior. These facts explain the raise of the effective "2a" energy 
with N, which means that the energy of the dilute Na states increases with and does not 
stay around the threshold energy. In Table HVI we can find quantitatively the increase of the 
effective 2a system with A^. The gradual dominance of the total Coulomb potential energy 
over the total 2a nuclear potential energy in the Na system can be also seen in Table HV1 

It is interesting to see the role of the 3a; potentials in the dilute Na system. From 
Table HVl the contribution to the total energy is small in the case of ^^C and ^^O, while it 
gives a non-negligible effect with increasing A^, as seen in Table ITTTl because it is proportional 
to the number of A^(A^ — 1)(A^ — 2)/6. In fact, the energy (V3) becomes steadily larger with 
the number of a bosons, although the quantity of the energy divided by the number of trios 
of the a bosons is decreasing with N, reflecting the fact that the A^a system is expanding 
with A^. It is instructive to study the case of no 3a potentials within the present framework. 
The calculated results show that the A^a system gradually falls into a collapsed state with 
increasing A^ as discussed before; for example, E' ~ — 55 MeV and {r%) = 4.77 fm for 
A^ = 40. The result tells us that the dilute a boson states appear under the cooperative 
work between the two-body and three-body potentials in the present model. 

Finally, we estimate the critical number N^r of the dilute A^a state beyond which the 
system is unbound. Since the critical state should be very unstable, it is difficult to determine 
it exactly with the present Hill- Wheeler-equation approach. However, we can deduce it 
approximately from studying the stability of the calculated eigen energies and rms radii 
against changing the model space determined by the Gaussian size parameters in Eqs. 
and The procedure is as follows. The present framework is able to describe dilute 

states trapped by the Coulomb barrier under the condition that we choose the range of the 
Gaussian size parameters wide enough to cover a configuration space over the whole rms 
radius of the states. If the state is stable or relatively stable, the calculated energies and 
rms radii are not changed very much against the variation of the Gaussian size parameters. 
Otherwise, those depend sensitively on the choice of the parameters. Using the wide model 
space, however, a special attention must be paid to investigate the eigen-states obtained 
from the Hill- Wheeler equation, because a discretized continuum state, which is unphysical, 
has a chance to become the lowest state in energy in the present variational calculation. We 
can easily identify it by investigating the behavior of its energy and/or rms radius against 
changing the Gaussian size parameters. The unphysical state usually has an abnormally 
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long rms radius, which is a httle smaller than or almost the same as the value of b\ """" or 
^{Kmax) (HH) and (123) • Table IVl shows the calculated results for the two parameter sets 

A and B (see Table ITTT|) . The calculated energies and rms radii for the dilute Na states with 
= 3 ~ 10 are almost the same for both the parameter sets A and B, while we see some 
discrepancies in those with A^ > 11, in particular, in the rms radii, for the two parameter 
sets. Even if we use other Gaussian size parameter sets which cover the configuration space 
over the rms radius for the dilute states, the results obtained are found to be similar to those 
mentioned above. Thus, the critical number of the dilute A^a boson systems is estimated 
roughly to be Ncr ~ 10, namely, ^°Ca, in the present study. It is noted that the number is 
the same as that in the spherical case as well as in the Gross-Pitaevskii-equation approach. 

IV. SUMMARY 

The dilute A^a cluster condensate states with J'^ = have been studied with the 
Gross-Pitaevskii equation (GP) and the Hill- Wheeler equation (HW), where the a cluster 
is treated as a structureless boson. Two kinds of effective a-a potentials were used: the 
density-dependent a-a potential (DD) and the phenomenological 2a potential plus 3-body 
a potential (PP), and we also included the folded Coulomb potential between the 2a bosons. 
Both potentials (only the latter) were (was) applied to study the structure of the spherical 
(spherical and axially deformed) dilute Na states with the Gross-Pitaevskii equation (Hill- 
Wheeler equation). Thus, we have studied the dilute multi a cluster states with the four 
different frameworks; (equation, potential, shape) = (GP, DD, S), (GP, PP, S), (HW, PP, S) 
and (HW, PP, D), where S and D denotes the spherical and deformed shapes, respectively. 

The main results to be emphasized here are as follows: 

1) All of the A^a states obtained show the dilute Na structure. They are common to 
all of the four different frameworks. The reason of why the total energy of the gas-like A^a 
state increases gradually with A^ and do not remain around their Na threshold values is 
understood as the competition between the nuclear a-a nuclear potential (attractive) and 
its Coulomb potential (repulsive). In fact, increasing the number of A^, the a-a distance is 
becoming larger (6 ~ 12 fm for A^ = 3 ~ 12), and then the contribution from the a-a nuclear 
potential per 2a pair in the dilute A^a state is decreasing rapidly because of the short-range 
character, while that from the Coulomb one is decreasing very slowly and remains almost 
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constant for = 5 ~ 12, reflecting tlie 1/r character of the Coulomb potential. 

2) The N dependence of the behavior of the calculated single a potentials obtained by 
solving the Gross-Pitaevskii equation within the (GS, DD, S) and (GW, PP, S) frameworks 
is impressive. For the 3a and 4a systems, we see an almost flat potential behavior in the 
inside region, and the Coulomb potential barrier in the outer region, while increasing the 
number of a bosons from = 5, the depth of the flat potential becomes shallower, and its 
range is expanding to the outer region, while the height of the Coulomb potential barrier 
is decreasing and almost disappears at around = 10. The origin of the appearance of 
the flat region is mainly due to the validity of the Thomas-Fermi approximation, namely, 
the neglect of the contribution from kinetic energy because of the dilute character of the 
A^a system. The result is analogous to the dilute atomic condensate state trapped by the 
magnetic flelds at very low temperature. On the other hand, the results for the total energy 
and rms radius were quite similar for the above-mentioned two frameworks except for A^=3 
and 4. The result might indicate that the dilute A^a states with A^ > 5 are not very sensitive 
to the details of the effective a-a potential (DD or PP). 

3) Comparing the results from the Gross-Pitaevskii equation (CP) and Hill-Wheeler 
equation (HW), we could see the effect of the neglect of the center-of-mass motion in the 
Gross-Pitaevskii equation and also the usefulness of the Hill- Wheeler-equation approach for 
describing the dilute A^a system. In the small-number A^a systems, the effect was found 
to be non- negligible, in particular, for the rms radius, while there is less than about 10 % 
deviation for the A^a > 5a systems. 

4) The axial deformation effect in the dilute A^a systems is substantial for the small- 
number a systems, but it is getting smaller for A^ > 6. In fact, the energy gain and 
reduction of the rms radius due to the deformation are 0.6 ~ 1.7 MeV and 11 ~ 19 %, 
respectively, for the former system. Due to the effect, the calculated results of the total 
energy and rms radius for A^=3 and 4 in the (HW, PP, D) framework are improved from 
those in the (HW, PP, S) framework, and give a good correspondence with those for the 
condensate a cluster states discussed by Tohsaki et al. llcJi] The present result that the dilute 
A^a states ( J'^ = 0+) may be deformed is natural, taking into account the fact that a gas-like 
A^a state with relatively small number of a's could easily be deformed. 

5) We estimated the critical number of the a bosons, N^r, beyond which the system is 
unbound for the four frameworks, (CP, DD, S), (CP, PP, S), (HW, PP, S) and (HW, PP, D). 
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All of the frameworks indicated that the number is roughly Ncr ~ 10, which is not strongly 
dependent on the frameworks. Thus, we concluded that the dilute Na cluster states could 
exist in the ^^C to ^°Ca systems with J'^ — 0+, whose energies vary from threshold up to 
about 20 MeV in the present calculation. 

The estimate of N^r ^ 10 is of course subject to the validity of our phenomenological 
approach treating the a-particles as ideal bosons. We, however, believe, for reasons outlined 
in the paper, that our estimate for Ncr may be correct to within ±20 %. In any case the 
value for Ncr, i-e. the maximum of a-particles in the condensate state is relatively modest. A 
very interesting question in this context is whether adding a few neutrons may stabilize the 
condensate and thus allow for much higher numbers of condensed a's. One should remember 
that *Be is (slightly) unbound whereas ^Be and ^°Bc arc bound. 

Concerning experimental detection of the a-condensates, the decay scheme of the dilute 
Na state is conjectured to proceed mainly via a decay. This indicates that such systems 
may be observed through the following sequential a decays: [dilute Na state] — > [dilute 
(N — l)a state]+Q;, [dilute {N — l)a state] [dilute {N — 2)a state]+Q;, • • •. Therefore, the 
sequential a decay measurement is expected to be one of the promising tools to search for 
the dilute multi-a cluster states, produced via the a inelastic reaction, heavy-ion collision 
reaction, and so on. It is highly hoped to perform such the experiments in near future. 
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TABLE I: Calculated results of the Gross-Pitaevskii equation with the density-dependent a- 
a potential in Eq. H29|): single a particle energy e and contributions from the kinetic energy 
(t), two-range-Gaussian term {v2g), density-dependent term (vd) and Coulomb potential (vc) in 
Eq. (|29j) . The total energy and nuclear rms radius for the Na system are denoted as E and 
respectively. All energies and rms radii are given, respectively, in units of MeV and fm. 



N 


nucleus 


e 


(t) 








E 


^) 


3 


12c 


0.18 


0.38 


-3.42 


2.31 


0.91 


0.38 


4.14 


4 


160 


0.78 


0.32 


-3.74 


3.04 


1.16 


1.44 


4.91 


5 


20Ne 


1.32 


0.28 


-3.90 


3.69 


1.25 


2.97 


5.53 


6 


24Mg 


1.82 


0.25 


-3.95 


4.25 


1.26 


4.94 


6.07 


7 




2.28 


0.22 


-3.91 


4.76 


1.21 


7.36 


6.58 


8 


32S 


2.72 


0.20 


-3.80 


5.20 


1.20 


10.2 


7.05 


9 


36 Ar 


3.13 


0.19 


-3.65 


5.59 


1.00 


13.4 


7.51 


10 


40Ca 


3.51 


0.18 


-3.45 


5.93 


0.86 


17.0 


7.98 


11 


44Ti 


3.87 


0.18 


-3.19 


6.19 


0.69 


21.0 


8.46 
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TABLE II: Calculated results of the Gross-Pitaevskii equation with the phenomenological la and 
3a potentials in Eqs. H31() and 1)32(1 : single a particle energy e and contributions from the kinetic 
energy (t), 2a potential (^2) in Eq. (|^T|) . 3a potential {v-^ in Eq. and Coulomb potential 
(fc)- The total energy and nuclear rms radius for the Na system are denoted as E and \J~{j^, 
respectively. All energies and rms radii are given, respectively, in units of MeV and fm. 



N 


nucleus 


e 


(t) 


{V2) 




^3) 


E 




3 


12c 


0.45 


0.25 


-1.90 


1.96 


0.14 


0.98 


4.87 


4 


160 


0.76 


0.27 


-2.71 


2.86 


0.33 


1.84 


5.23 


5 


20Ne 


1.12 


0.27 


-3.36 


3.68 


0.53 


3.04 


5.55 


6 


24Mg 


1.51 


0.26 


-3.89 


4.44 


0.70 


4.63 


5.85 


7 


28Si 


1.91 


0.26 


-4.31 


5.13 


0.84 


6.61 


6.13 


8 


32s 


2.31 


0.24 


-4.64 


5.78 


0.94 


8.99 


6.40 


9 


36 Ar 


2.73 


0.23 


-4.89 


5.37 


1.02 


11.8 


6.68 


10 


40Ca 


3.13 


0.22 


-5.06 


6.91 


1.60 


15.0 


6.95 


11 


44rpi 


3.53 


0.20 


-5.15 


7.40 


1.70 


18.6 


7.24 
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TABLE III: Gaussian size parameter sets A and B used to obtain the deformed Na states with 
the axial symmetry by solving the Hill- Wheeler equation. The unit of h\ and 63 is fm. 



set 


"'max 


6?) 






"3 




A 


12 


0.4 


9.0 


12 


0.5 


9.5 


B 


12 


0.4 


11.0 


12 


0.5 


10.5 
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TABLE IV: Calculated total energy E, total kinetic energy (T), total 2a nuclear potential energy 
(V2), total Coulomb potential energy {Vc), total 3a potential energy (V3) and nuclear rms radius 
^ (r^) for each dilute Na state. The effective energy of the "2a" system in the A''a system is 
denoted as £2 (see text). All energies and rms radii are given in units of MeV and fm, respectively. 









deformed 


case 








spherica' 


case 




nucleus 


(T) 




m 


(^3) 


E 


(^2) 






/K) 


3 


12c 


2.24 


-5.95 


3.62 


0.09 


-0.01 


(-0.00) 


3.73 


0.64 


4.18 


4 


160 


3.74 


-11.31 


7.21 


0.42 


0.11 


(0.02) 


3.90 


1.58 


4.74 


5 


20Ne 


4.32 


-15.92 


11.58 


1.13 


1.11 


(0.11) 


4.20 


2.83 


5.19 


6 


24Mg 


3.71 


-18.78 


16.05 


2.09 


3.13 


(0.21) 


4.69 


4.47 


5.64 


7 


28s 


3.13 


-21.03 


20.58 


0.09 


5.58 


(0.27) 


5.24 


6.48 


5.99 


8 


32Si 


2.77 


-23.16 


25.31 


0.09 


8.30 


(0.30) 


5.79 


8.93 


6.33 


9 


36 Ar 


2.52 


-25.11 


30.22 


0.09 


11.31 


(0.31) 


6.34 


11.81 


6.90 


10 


40Ca 


2.35 


-26.76 


35.21 


0.09 


14.62 


(0.32) 


6.90 


14.98 


7.26 


11 




2.26 


-28.52 


40.58 


0.09 


18.27 


(0.33) 


7.38 


18.53 


7.54 


12 


48Cr 


2.21 


-30.43 


46.37 


0.09 


22.63 


(0.34) 


7.78 


22.44 


8.62 
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TABLE V: Calculated energies E of dilute Na states (J^ = 0"*") together with the nuclear rms 
radii (r^) and rms distances between 2a bosons \/ ('^ao), where we use the two Gaussian size 
parameter sets A and B, shown in Table ITTTI The energy E is measured from the respective Na 
threshold. The units of energy and rms radius (distance) are MeV and fm, respectively. 



N 


nucleus 


E 


set A 




E 


set B 




3 


12c 


-0.01 


3.73 


5.73 


-0.01 


3.72 


5.72 


4 


160 


0.11 


3.90 


5.72 


0.11 


3.90 


5.72 


5 


20Ne 


1.11 


4.20 


6.06 


1.13 


4.20 


6.06 


6 


24Mg 


3.13 


4.69 


6.76 


3.13 


4.68 


6.75 


7 


28Si 


5.58 


5.24 


7.57 


5.59 


5.24 


7.57 


8 


32s 


8.30 


5.79 


8.36 


8.30 


5.78 


8.35 


9 


36 Ar 


11.31 


6.34 


9.16 


11.31 


6.35 


9.17 


10 


40Ca 


14.62 


6.90 


9.97 


14.61 


7.03 


10.58 


11 




18.27 


7.38 


10.65 


18.22 


7.89 


11.42 


12 


48 


22.27 


7.78 


11.21 


22.10 


8.81 


12.77 
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FIG. 1: Radial part of the relative wave function between the 2a clusters in the = 0^ resonant 
state at i?2a=92 keV with use of the soft-core q-q potential in Eq. pif) . The scale of r(j){r) is 
arbitrary. 
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FIG. 2: (a) Total energies for the dilute Na states measured from each Na threshold, and 
(b) their nuclear rms radii, which are obtained by solving the Gross-Pitaevskii equation with the 
density-dependent potential. 
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FIG. 3: Single a particle potentials Ua{R) (solid line) which are obtained by solving the Gross- 
Pitaevskii equation with the density-dependent potential; (a) 3a, (b) Aa, (c) 5a, (d) 6a, (e) 7a, 
(f) 8a, (g) Qa, (h) lOa, (i) 11a systems. The dashed, dot-dashed and dotted lines demonstrate, 
respectively, the contribution from the two-range Gaussian term, density-dependent term and 
Coulomb potential. 
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FIG. 4: (a) Total energies for the dilute Na states measured from each Na threshold, and 
(b) their nuclear rms radii, which are obtained by solving the Gross-Pitaevskii equation with the 
phenomenological 2a and 3a potentials. 
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FIG. 5: Single a particle potentials Ua{R) (solid line) which are obtained by solving the Gross- 
Pitaevskii equation with the phenomenological 2a and So; potentials; (a) Sa, (b) 4a, (c) 5a, 
(d) 6a, (e) 7a, (f) 8a, (g) 9a, (h) 10a, (i) 11a systems. The dashed, dot-dashed and dotted 

lines demonstrate, respectively, the contribution from the 2a potential, 3a potential and Coulomb 
potential. 
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FIG. 6: (a) Total energies for the dilute Na states (spherical case) measured from each Na 
threshold (solid line), and (b) their nuclear rms radii (solid), which are obtained by solving the 
Hill- Wheeler equation with the phenomenological 2a and 3a potentials. For comparison, we also 
give the corresponding values obtained by solving the Gross-Pitaevskii equation with use of the 
same potentials (dotted lines). 
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FIG. 7: (a) Total energies for the dilute Na states with the axial deformation measured from each 
Na threshold (solid line), and (b) their nuclear rms radii (solid) and rms distance between 2a par- 
ticles (dotted) , which are obtained by solving the Hill- Wheeler equation with the phenomenological 
2a and 3a potentials. 
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FIG. 8: Phenomenological 2a potential (solid line), where the 2a nuclear potential [Eq. pij) ] and 
its Coulomb potential [Eq. pUj) ] are drawn by the dashed and dotted lines, respectively. 
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